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We extend an earlier study of the Landau gauge quark propagator in quenched QCD where we 
used two forms of the C(a)-improved propagator with the Sheikholeslami-Wohlert quark action. 
In the present study we use the nonperturbative value for the clover coefficient Csw and mean- 
field improvement coefficients in our improved quark propagators. We compare this to our earlier 
results which used the mean-field Csw and tree-level improvement coefficients for the propagator. 
We also compare three different implementations of tree-level correction: additive, multiplicative, 
and hybrid. We show that the hybrid approach is the most robust and reliable and can successfully 
deal even with strong ultraviolet behavior and zero-crossing of the lattice tree-level expression. 
We find good agreement between our improved quark propagators when using the appropriate 
nonperturbative improvement coefficients and hybrid tree-level correction. We also present a simple 
extrapolation of the quark mass function to the chiral limit. 
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I. INTRODUCTION 



Lattice studies of the quark propagator provide a direct, model-independent window into the mechanism of dynam- 
ical chiral symmetry breaking and its momentum dependence. In addition it provides insight into the the nature and 
location of the transition region of QCD where inherently nonperturbative behavior evolves into the more analytically 
accessible perturbative form. Furthermore, direct lattice calculations of the quark propagator inform hadron model 
building from the relativized constituent quark picture to quark models based on Schwinger-Dyson equations j^, |^] . 

In a recent paper |^ we presented a method for removing the dominant ultraviolet tree-level lattice artifacts in 
the momentum-space quark propagator. This was a generalization of the concept of tree-level correction, which 
was first introduced in the study of the gluon propagator |4[ |[ ^, 0]. It was shown that, for two different 0{a)- 
improved propagators Si and Sr, see Eqs. (||) and and using a mean-field improved action, this leads to a 
dramatic improvement in the ultraviolet behavior of the propagator. However, the remaining ultraviolet artifacts are 
sufficiently large to make the results unreliable beyond pa ~ 1.2. Moreover, the two improved propagators remain 
discernibly different even in the infrared, yielding different estimates for the infrared quark mass. 

Here we will present results using nonperturbatively determined values for the 0{a) improvement coefficients, rather 
than the tree- level and mean- field improved coefficients used in Ref . |^ . We will also present two alternative techniques 
for removing tree-level artifacts and will discuss the relative merits of the three methods. 



II. IMPROVEMENT 



The general scheme for 0{a) improvement of the quark propagator was discussed in Ref. Here we restrict 
ourselves to presenting the formulae and definitions which we will be using in this paper. For further details, see 
Ref. [Pl and references therein. 
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The SW fermion action, 

£{x) = £^{x) - ^Cswa'ipix)(Tf_,„Ff_,„{x)il;{x) , (1) 

combined with appropriate improvements of operators can be shown |^ to remove all 0{a) errors in on-shell matrix 
elements. For ofF-shell quantities such as the quark propagator it is not that simple, and no general proof of 0{a) 
improvement is known. Indeed, to calculate gauge dependent quantities one might expect to have to introduce gauge 
non-invariant (but BRST invariant) terms in the action. However, at tree level it is possible to proceed by adding all 
possible dimension-5 operators to the action and eliminating all but the clover (SW) term by a field redefinition 
Beyond tree level, one may proceed by adding all possible terms with the correct dimensionality and quantum numbers 
to the operator in question, and tuning the parameters to eliminate 0{a) terms. Ignoring the gauge non-invariant 
terms (which are discussed in Refs. |ll|) we may write down the following expressions for the 0{a) improved 
quark propagator, 

Si{x, y) = {Si{x, y; U)) = ((1 + bqam)SQ{x, y; U) - aX5{x - y)) , (2) 
SRix,y) EE {SR{x,y;U)) = {{l + b'^am) [l - c'^(x)] So{x,y;U) [l + c'^f (y)]) . (3) 

where So{x, y; U) for a given configuration U is the inverse of the fermion matrix. Note that if we are only interested in 
on-shell quantities such as hadronic matrix elements, the (5-function can be ignored, so we only need 5o together with 
the improvement coefficients for the various operators. From this, 5/ is easily obtained, whereas Sr is computationally 
somewhat more expensive. 

The coefficients hq, b'^, A and c'^ must be tuned in order to eliminate 0{a) errors in the propagator (as far as this 
is possible). At tree level, their values are &g = 1, A = 6^ = i, = i. The values for hq and A have recently been 
calculated at one-loop level |Q . The mean-field improved values for all these coefficients may be obtained by dividing 
the tree-level values by the mean link uq. The one-loop mean-field improved values have also been calculated in 
Ref. but as we will argue, the small changes in values obtained by including the one-loop contribution have very 
little practical effect, so we will not use these here. It should also be noted that the mean-field improved value for A 
is very close to the nonperturbative value reported in Ref. . This indicates that mean-field improvement of the 
coefficients bq and A (or, alternatively, h'q and c^) may be sufficient to remove 0{a) errors to the desired precision. 

The bare mass also receives an 0{a) correction, which can be expressed as follows 

m = (1 + b„iam)m , am = . (4) 

The coefficient has been calculated at one- loop order [Q, 

5„, = -i- 0.09625g + O(.g4). (5) 
When evaluating Eq. (^, we will be using the boosted coupling constant = 9o/'^o- 

III. TREE-LEVEL CORRECTION 

In the continuum, the spin and Lorentz structure of the quark propagator, together with parity symmetry, deter- 
mines that the propagator must have the following form, 

c( \ ^ Z{fi;p'^) ^ 1 

On the lattice, what we measure is the bare (regularized but unrenormalized) propagator. This differs from the 
renormalized propagator in Eq. (|6|) by an overall renormalization constant Z2(/^,a), which we will absorb into Z{p), 
as we did in Ref. to simplify the presentation of our results. 

In Ref. ^ we defined a tree-level correction procedure involving an overall multiplicative correction and an additive 
correction of the mass function, as follows, 

S-\pa) = z{pa)Zi^){pa) ^ ^ + ^M^iv^) + aAM(°)(pa)] , (7) 
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where fc^ = sin(p^a)/a and Z^"^ and AM^^^ are defined by the tree-level quark propagator, 

i + am + aAM^°'^ {pa) . (8) 



(s^^Hpa))' 



Z^^){pa) 



The functions Z(j)a) and M"'{pa) should then express the nonperturbative behavior of the quark propagator, with 
the dominant lattice artifacts removed. We saw that this procedure led to a dramatic improvement in the behavior 
of these functions, but at large momenta the data could still not be trusted because of large cancellations. 
Here we will consider an alternative, purely multiplicative tree-level correction procedure, defined by 

^"'^^"^ " Z{pa)m{pa) h ^ + , (9) 

where Zm'' (pa) is defined by the tree-level expression, 

(5(°)(pa))"' = ^^j^ [z ^a + amZ^^\pa)\ . (10) 

The tree- level corrected mass function M is thus obtained from the uncorrected function = tr /ANc via 

aM™(pa) = M^(pa)/Z(?)(pa) . (11) 

This procedure should not suffer from the problem of large cancellations. However, it will encounter problems when 
either the numerator or denominator of Eq. (^l|) crosses or is close to zero. In order to remedy this problem, we consider 
a third, 'hybrid' scheme, where the negative part of the tree-level expression is subtracted, while the remaining positive 
part is multiplicatively corrected. Specifically, we define AM'+', AM'"' such that 

AM'+' (pa) + AAf (pa) = AM^") (pa) (12) 

AM<+'(pa)>0; AM<-'(pa)<0 Vpa (13) 

Then we can write 

M(") (pa) ^ am + AM<+> (pa) -I- AAf <-> (pa) = amZ<+' (pa) + AAf <-> (pa) . (14) 

The tree- level corrected mass function M^{pa) is then 

aM^{pa) = {M^ipa) - aAA'p-\pa)) lZ^^\pa) . (15) 

The definition of this scheme contains an ambiguity, since it is obvious that we may still satisfy (|l^), ( p^ by adding 
any strictly positive term to AAf*"*"' and subtracting the same term from AAf*~'; e.g. by taking AAf*+' AAf'+' + 
AAf AAf'-' - k^. In order to remove this ambiguity, we add the following criteria: 

1. Factoring out a common (positive) denominator, AAf '+' and AAf '"' should be polynomials in the 4 variables 
A:^, A:^, Afc^ and m, where we have defined 

2 

/fc^ = - sin(p^a/2); a^Afc^ ee fc^ - /fc^ . ^;^g) 

2. The coefficient of each term must be positive for AAf<+' and negative for AAf'"'; 

3. Any one monomial in fc^,fc^,Afc^ and m can only occur in one of AAf<+' or AAf'"'; eg., if there is a term 
proportional to mfc^ in AAf <"' there cannot be a term proportional to mfc^ in AAf 

These criteria ensure that AAf <+' and AAf'"' are as small as possible, leading to the minimum possible distortion of 
the data. 

Specifically, the expressions we use are 

, , , , , am (&„ - A)a2m2 + Aa^'AP + (i& A)a3mfc2 , , 

^™ /Op = — — 7T- — ' 

' am(l -I- am) 

„AMr(ri^-i^*4±|^**, (18) 

2(1 + am) 



aAAf'"" 
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Volume 


a ^ (GeV) 


Csw 


K 
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am 


m (MeV) m 


(MeV) iVcfg 


6.0 


16^ X 48 


2.120 


MF 


0.13700 


Si 


0.0579 


123 


118 499 












Sr 






20 










0.13810 


Si 


0.0289 


61 


60 499 


6.0 


le^* X 48 


2.120 


NP 


0.13344 


Si 


0.0498 


105 


102 10 












Sr 






20 










0.13417 


Si 


0.0294 


62 


61 10 










0.13455 


Si 


0.0188 


40 


39 10 


6.2 


24'' X 48 


2.907 


MF 


0.13640 


Si 


0.0399 


116 


113 54 










0.13710 


Si 


0.0212 


62 


61 54 



TABLE I: Simulation parameters. 'MF' and 'NP' refer to the mean-field improved and nonperturbatively determined values 
for Csw respectively. The improved propagators Si and Sr are defined in Eqs. (S) and (0). 



where we have written 

AUp) = 1 + -am + —a'^k^ + -a'^AP . (21) 
2 16 4 

(22) 

It should be remembered that the mass function M{p) must be renormalization-point independent in a renormal- 
izable theory and that the current quark mass at the renormaUzation point m{fi) is given by m(/i) — M{p = fi). 
The ultraviolet mass function is of course only constant up to logarithmic corrections. The multiplicative and hy- 
brid tree-level correction ensures that the zeroth-order perturbative behavior of the mass function in the ultraviolet 
matches that of the continuum. The logarithmic corrections should principle show up in the lattice data, as they did 
for the gluon propagator in Ref. However, this is a small effect compared to the tree-level lattice artifacts. It will 
be a measure of the success of our improvement and correction scheme whether the logarithmic corrections may be 
extracted from the lattice data. 



IV. RESULTS 

In addition to the data used in Ref. we have analyzed data at (3 — 6.0 using the nonperturbatively determined 
value for Csw (=1.769), and at /3 = 6.2 using the mean-field improved Csw (=1.442). This is a subset of the UKQCD 



data analyzed in Ref. |15 , with the values for Csw taken from Ref. pi| . The simulation parameters are given in 
Table |. Note that the values of am are different from those given in Ref. ||^; this is because we have here used 
the determination of Hc reported in Ref. p5[ | instead of an earlier, preliminary value. All the data shown have been 
obtained from the raw data using the cylinder cut described in Ref. The scale is taken from the hadronic radius 
To [0 using the interpolating formula of Ref. ||l^ and the phenomenological value tq = 0.5 fm. Note that this differs 
from the scale used in Ref. Jp| , which was taken from an earlier determination of the string tension. The gauge fixing 
is identical to that of Ref. |^ which was also used in Ref. ||^ . This is a version of lattice Landau gauge that contains 
Gribov copies; the effects of these have not been studied here. 



A. Results with mean-field improved Cs 



We first consider the effect of employing the multiplicative and hybrid correction schemes on the data analyzed in 
Ref. Figure Q shows the tree-level corrected mass function M evaluated using the three schemes, for both Si and 
Sr at f3 — 6.0, Csw = MF. It is clear that the ultraviolet behavior is much improved, but the multiplicatively corrected 
M from Si exhibits pathological behavior at intermediate momenta. This is a consequence of a zero crossing in the 
tree-level mass function, leading to division by near-zero numbers in Eq. (pT|). In the hybrid scheme, this problem is 
absent. 

It is worth noting that although the mass function approaches the subtracted bare mass m in the ultraviolet, the 
actual values obtained using the multiplicative and hybrid schemes differ from each other and from the bare mass by 
up to 20%. It is clear that at this stage this procedure is not good enough to yield a good estimate of current quark 
masses. We also see that there is no sign in these data of the logarithmic running of the current quark mass. 

For Sr, we also see a clear improvement in the ultraviolet behavior, as well as a small but significant difference in 
the ultraviolet mass between the multiplicative and hybrid schemes. Since the tree-level mass function for Sn does 
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FIG. 1: The tree-level corrected mass function, for Si (left) and Sr (right) with Csw ~ MF, k — 0.13700, using additive, 
multiplicative and hybrid correction. The hybrid scheme is robust even in the presence of a zero crossing of the tree-level mass 
function for Si. 



not have any zero crossings, the multiplicatively corrected mass does not exhibit the same pathological behavior as 
for Sj. The tree-level mass function for Sr does not cross zero, but does approach it, and the effects of this may be 
detected at intermediate momenta. Hence we consider the hybrid correction to be more reliable for Sr. 



B. Results with nonperturbative c. 



We now turn to the effect of using nonperturbatively determined, rather than tree-level or mean-field, improvement 
coefficients. Figure ^ shows, on the left, the mass function obtained from our two improved propagators using the 
nonperturbative value for Cgw and the mean- field values for bq^b'^ and A (we have not been in a position to obtain 
data using mean-field or nonperturbative values for c^, only the tree- level value c'^ = j). As indicated in Sec. II, the 
nonperturbative values for the latter coefficients are currently not known, but at least the nonperturbative value for 
A reported in Ref. Q is close to the mean-field value, and it seems reasonable to guess that this is the case for the 
other coefficients as well. We therefore assume that although it is not entirely consistent, it is not unreasonable to use 
the mean-field values. On the right of Fig. || are shown equivalent data from Ref. |^ , along with data using mean-field 
values for all the improvement coefhcients in Sj. 

It is immediately clear that using the NP value for Csw removes the large discrepancy between the two improved 
propagators for pa < 2, even when using additive tree-level correction. For larger pa, the discrepancy remains, which 
is not unexpected since at these momenta 0{a^) and higher errors become dominant. 

It is instructive to compare this with the effect of reducing the lattice spacing, as shown in Fig. ^. We do not have 
data for Sr at P — 6.2, but we see that Si changes very little with /3 in the intermediate momentum range where 
the discrepancy becomes large. Assuming that Sr, changes by a similar amount, we may conclude that reducing the 
lattice spacing does very little to reduce the discrepancy, although the behavior of Si at large momenta is somewhat 
improved, as one would expect. We can also see that going from tree-level to mean-field values for the coefficients bq 
and A has only a very small effect. We may also conclude that using one-loop values for these coefficients will have 
negligible effect, since the difference between tree-level and one-loop coefficients is even smaller than that between the 
tree-level and mean-field improved values. This also gives us added confidence in the use of mean-field rather than 
nonperturbative values for these coefficients. 

In Fig. ^ we show the tree-level corrected Z{p) function, for both nonperturbative and mean-field Csw The upper 
figures show the 'unrenormalized' Z{p) — the upper right figure is taken directly from Fig. 5 in Ref. We see 
that there is still a very significant discrepancy between Si and Sr, even with the nonperturbative Cgw Much of this 
discrepancy, however, amounts to an overall renormalization, which may be included in the quark field renormalization 
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FIG. 2; The additively tree-level corrected mass function, for c^w = NP, k = 0.13344 (left) and c^w = MF, k = 0.13700 (right 
S/,imp in the right-hand figure is obtained using the mean-field improved rather than the tree-level values for bq and A in Eq, 
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FIG. 3: The additively tree-level corrected mass function, for Csw = MF, at /J = 6.2 and f3 = 6.0. 



constant Z2. To eliminate this possible, unphysical source of disagreement we rescale the data by imposing the 
'renormalization condition' Z{pa = 1) = 1. The result of this is shown in the lower panel of Fig. ^. We then see that 
the infrared behavior of the tree-level corrected Z{p) functions agree much better than they did in the previous work 
in Ref. For the few most infrared points we see that there is an apparently better agreement for Z{p) between the 
two forms of the propagator when the nonperturbative Cgw is used. The ultraviolet agreement is less than satisfactory 
even when the nonperturbative Csw is used and we conclude that the ultraviolet behavior of the tree-level form of Si 
is too severe to be remedied by our tree-level correction scheme even with nonperturbative and mean-field improved 
coefficients in the action and propagators respectively. Because of the more reasonable tree-level behavior of Sr and 
because in this case Z{p) is almost unchanged when using either the nonperturbative or mean- field Cgw, we take as 
our best estimate for Z{p) the tree- level corrected result from Sr with nonperturbative Csw 
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FIG. 4: The tree-level corrected Z(p), for c^„ = NP, k = 0.137 (left) and Csw = MF, k = 0.13344 (right). The upper two 
figures show the data in lattice units, and without rescaling. The lower two figures show Z{p) vs. momentum in physical units, 
and after rescaling ( "renormalizing" ) so that ^(2.1GeV) — 1. The infrared agreement after rescaling is very good and we take 
the tree-level corrected Z{p) from Sr with the nonperturbative Csw as the best estimate for this quantity (see text). 



Figure || shows the mass function for Csw = NP, with the multiplicative and hybrid correction schemes. The 
multiplicative scheme exhibits the same problems as those we encountered with Cg^^ = MF. In particular, Zm"^ (p) 
from Si crosses zero for pa ^ 0.5 and renders the multiplicative scheme meaningless for that case. Even for Sr we 
find that the uncorrected has small zero-crossings for momentum values in the range 1.5 < pa < 2.3, which 
render the multiplicative scheme unsatisfactory. Using the hybrid scheme, however, we avoid these pathologies of the 
naive multiplicative scheme as can be seen from the two lower figures in Fig. ||. In Fig. ^ we have plotted the mass 
functions for the two propagators in physical units. We see good agreement between the mass functions using hybrid 
tree-level correction and the nonperturbative Cgw coefficient across the entire range of available momenta. The mass 
function for Sr dips slightly below that for Sj at intermediate momentum points even though they approach very 
similar asymptotic values. This residual disagreement implies that we have not succeeded in removing all of the lattice 
artifacts at intermediate momenta, although hopefully we have gone some significant way toward achieving that end. 
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FIG. 5: The tree-level corrected M{p), for Csw = NP, using 5*/ (left) and Sr (right), and with the multiplicative (top) and 
hybrid (bottom) correction schemes. The multiplicative scheme clearly fails for 5*/, and also performs poorly for Sr; while the 
hybrid scheme performs well for both and leads to good agreement between the two propagators. 



C. Chiral extrapolations 



We have available data for three quark masses for Sj with the nonperturbative Csw Having seen the very plausible 
behavior of the mass function for 5/ after hybrid tree- level correction and the good agreement with that for Sr, we 
are given the confidence to attempt a simple extrapolation of the quark mass function to the chiral limit. The first 
chiral extrapolation we performed was a linear extrapolation of the ultraviolet mass (obtained by fitting M^{p) to a 
constant in the range 2 < pa < 3) as a function of fh. This is shown in Fig. where we see that the ultraviolet mass 
vanishes in the chiral limit to a very good approximation as it should. Indeed, the extrapolated value of —1 MeV 
is much smaller than the systematic uncertainties arising from the different tree-level correction schemes discussed 



above. We also see that for the ultraviolet mass a linear extrapolation does very well. As we noted in Sec. Ill, the 
ultraviolet mass function is only constant up to logarithmic corrections, but our ultraviolet behavior is not sufficiently 
under control that it would be meaningful to attempt to extract those from Fig. ^. 

Also in Fig. |^ we show the result of a linear extrapolation of the infrared quark mass M{p — > 0), together with our 
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FIG. 6: The hybrid tree-level corrected M^{p), for c^w = NP, using Si (left) and 5"^ (right). For m = 91 MeV we find good 
agreement between the two data sets, both in the infrared and ultraviolet. The residual disagreement at intermediate momenta 
is a pointer to lattice artifacts that we have not brought under full control, even with nonperturbative improvement and hybrid 
tree-level correction. 



previous results from Ref. Due to the small statistics, the error bars for the nonperturbative Cs„ data are quite 
large. We still observe that the extrapolated mass value for the nonperturbative Csw is systematically lower than for 
the mean- field Csw, but the results for Si and Sr are now fully consistent, and also agree with the value obtained 
from Sr with the mean- field Csw Note that the value in MeV for M{p — > 0) from the mean- field Csw differs from that 
reported in Ref. This is due to the different values for the lattice spacing in the two papers. The uncertainty in 
the lattice spacing adds an additional uncertainty of about 10% to all numbers in physical units. This uncertainty is 
an intrinsic feature of the quenched approximation. 

We also show in Fig. ^ the result of a simple quadratic chiral extrapolation in m for the entire mass function (for 
Si with hybrid tree-level correction and nonperturbative Csw)- The result has a plausible form and is presumably 
a good estimate of M{p) in the chiral limit. The very small dip at p ~ 1.4 GeV is within two standard deviations 
and, while not statistically significant, it is a again a hint that we have not completely removed lattice artifacts at 
intermediate momenta. A linear chiral extrapolation, while adequate in the ultraviolet and infrared, does not fit the 
data in the intermediate momentum regime. We see that the quadratic extrapolation to the chiral limit is consistent 
with a vanishing current quark mass and a rapid falling off of M{p) such that it appears to essentially vanish by 
approximately 1.5 GeV. This suggests that the effects of dynamical chiral symmetry breaking become negligible at a 
scale which we estimate to be = 1.45ii3 GeV, where the errors are purely statistical. 

We have also studied the systematic uncertainties arising from the specific choice of tree-level correction scheme. 



In order to quantify this, we have modified the hybrid scheme defined in Sec. [11 by taking 



AAf<+>(pa) ^ AA/<+>(pa) + eAA/<-'(pa) AAf<-'(pa) ^ (1 - e)AAf<-'(pa) (23) 

where e is a free parameter. Considering small variations in e, —0.1 < e < 0.1, we find that the correction scheme 
dependence gives rise to uncertainties in p-^ of about 100 MeV. Using the mean-field Csw at /3 = 6.0 and 6.2, we get 
values for p^ that are slightly higher, but still consistent withint two standard deviations. 

We take as our best estimate for the chiral symmetry scale the value from 5*/ with the nonperturbative Csw at /3 = 6.0: 
Px = 1.451*11318(14) GeV, where the first set of errors are statistical, the second are the systematic uncertainties due to 
the tree-level correction scheme, and the third is the uncertainty in the lattice spacing. This value is roughly consistent 
with the chiral symmetry breaking scale A^^sb arising in low-energy effective theories and in instanton models (see 
e.g. Re f. |19| ). An understanding of the relationship between this result and the recent analysis of the pseudoscalar 
vertex pO| is an interesting topic for future investigation. 
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FIG. 7: Left: The ultraviolet quark mass Muv as a function of the bare quark mass m, for Si with Csw = NP and using hybrid 
tree-level correction. The values are obtained by fitting M^{p) to a constant for 2 < po < 3. Right: The infraxed quark mass 
Mir = M{p = 0), obtained by extrapolating M{p) to pa = 0, as a function of rh. The bursts indicate the chirally extrapolated 
values of Mir obtained by a simple straight line fit for each action. The solid line represents the fit for Csw = NP, while the 
dotted line is the fit for Csw = MF. We see that the values of Mir from Si and 5^ agree for Csw = NP, giving us further 
indication of the superiority of nonperturbative to mean-field improvement, despite the large statistical uncertainties in the NP 
data. 
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FIG. 8: The hybrid-corrected mass function from Si with Csw = NP, with the bare mass rh extrapolated to zero using a 
quadratic fit. The small dip at p ~ 1.6 GeV is not statistically significant and may be due to residual lattice artifacts. The 
non-zero values for M(p) in the chiral limit are entirely due to dynamical chiral symmetry breaking and provide a direct measure 
of this effect. 
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V. DISCUSSION 



We have made use of asymptotic freedom to factor out the dominant (tree-level) lattice artifacts in the quark 
propagator at high momenta. We have discussed several different schemes for applying this idea, referred to as tree- 
level correction, to the mass function, which is the scalar part of the inverse quark propagator. A purely multiplicative 
scheme is seen to encounter problems because the value of the tree-level mass function approaches or crosses zero, 
leading to ill-defined behavior at intermediate momenta. The purely additive scheme defined in Rcf. although 
leading to a dramatic improvement on the uncorrected data, did not give reliable results for the mass function above 
pa ^ 1. We have defined a hybrid tree-level correction scheme which combines the additive and multiplicative schemes 
in such a way that the mass function becomes well-behaved at all momentum values. Ambiguities in the correction 
scheme should show up most clearly at intermediate momenta. We find that the uncertainties due to such ambiguities 
are of comparable magnitude to the statistical uncertainties. 

We have also studied the effect of using the nonperturbative value for Cs„. This was seen to improve the data 
considerably, as demonstrated most dramatically by the very good agreement between the two definitions Si and Sr 
of the improved quark propagator. In contrast, reducing the lattice spacing by going from P — 6.0 to /3 = 6.2 with a 
mean- field Cgw only gave a slight improvement. 

Finite volume effects may be estimated by studying the spread of points with momenta in different directions in 
the infrared, as was done in Refs. |^. We do not find any significant anisotropy at low momenta either in Z(p) or 
in M{p), so we conclude that finite volume effects here are small. 

We have not here considered the possible effect of Gribov copies. This remains an interesting subject for future 
study, which is currently being pursued. 

In summary, the key results of this study are that one should use the most appropriate (nonperturbative) determi- 
nations of improvement coefficients wherever possible. Where tree-level behavior is severe with zero-crossings or near 
zero-crossings the hybrid tree-level correction scheme can be used in place of the multiplicative one. In the infrared 
and intermediate momentum regimes we appear to have reasonable control over lattice artifacts. The chiral behavior 
resulting from a simple chiral extrapolation appears reasonable. We believe that the best estimate of the continuum 
Z{p) function corresponds to the nonperturbative Cgw result for the Sr propagator as shown in Fig. |[ All of the 
hybrid corrected mass functions with nonperturbative Csw appear reasonable and agree with each other. The chiral 
extrapolation of these is shown in Fig. ^. We find that the effects of dynamical chiral symmetry breaking become 
negligible above a momentum scale p^. Our best estimate for this scale is p-^ = 1.4511311(14) GeV. 

We emphasize that the real test of these conclusions will be to implement these methods on finer lattice spacings, 
with further improved actions, and ideally with actions which respect chiral symmetry on the lattice. These subjects 
are being pursued. 
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